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ABSTRACT 

The transition of a standard thin disk to a radiatively inefficient accretion flow (RIAF) is 
expected to occur, when m ~ m C rit (m = M/MEdd)- The radiative efficiencies of accretion 
flows accreting at rates lower than the critical accretion rate m cr i t become significantly lower 
than that of standard thin disks. It is believed that the initial transition radius is small just 
after the accretion mode transition, and then the transition radius increases with decreasing 
accretion rate, as suggested by some theoretical models and observations. Based on such 
variable transition radius models, we derive how the accretion rate rh(t) evolves with time 
from the observed Eddington ratio distribution for a sample of low-luminosity active galactic 
nuclei in the local universe. The derived time-dependent accretion rates rh(t) show a rapid 
decrease after the transition of the standard thin accretion disk to a RIAF, which is consistent 
with that derived from the hard X-ray background. 
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1 INTRODUCTION 

Accretion onto massive black holes is believed to power active 
galactic nuclei (AGNs). The UV/optical continuum emission ob- 
served in luminous quasars is attributed to the thermal radiation 
from the acc retion disks surroun ding the massive black holes in 
quasars (e.g., [Sun & Malkan 1989). In recent years, some different 
approaches are propos ed to measure th e masses of the central black 
holes in AGNs (e.g ., iPetersorj 1 19931 : iFerrarese & Merrittl l200d : 
iGebhardt et alJ|2000T) . and the central black hole masses of many 
AGNs can be measured fairly accurately. It is found that a frac- 
tion of luminous AGNs are accreting at extremely high rates. Their 
bolometric luminosities are around (or even higher than) the Ed- 
dington luminosity, for example, the black holes in many narrow- 
line Seyfert 1 galaxies are belie ved to be accreting at Eddington 
(or super-Eddington) rates (e.g.JSulentic et aL|l2000l ; IWarner et al.l 
l2004l : bian & Zhadl2004k IChen & Wandl2004l) . Besides the lumi- 
nous AGNs accreting at high Eddington rates (> 0.01), many 
nearby low-luminosity AGNs are found to have similar characteris- 
tics as those luminous AGNs, but with relatively weaker b road-line 
emiss ion and lower luminosities in different wavebands dHo et al.l 
1997). Most of these low-luminosity AGNs are accreting at highly 
sub-Eddington rates (Lboi/iEdd <C 0.01). The Eddington ra- 
tios of AGNs can spread over more than ten orders of magni- 
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tude (from < 10~ 10 to more than unity)(e.g., IWu & Cao|[2005l : 
IHopkins. Naravan & Hernquistll2006|) . 

It is believed that standard thin disks (or slim disks) 
are present in luminous AGNs, while the radiatively in- 
efficient accretion flows (RIAFs) are suggested to be 
present in those low-luminosi t y AGNs accreting at rela - 
tively low rates l Naravan & Yil 19941: Blandford & Begelma 



ely 



1999, 2004; Igumenshchev, Abramowicz & Narayan 



Naravan. Ig umenshche v & AbramowicauOOOl) . There is a critical 
accretion rate m cl -it, above which the RIAF is suppressed and 
a standard thin disk is present. The RIAF may connect to a 
standard thin disk at a certain transition radius R tI . This is 
required by modelling on a variety of ob s ervations of AGNs (e.g., 



Ouataert et all 1 19991; iLu & Wand l200d : Icaol 120031 ; iKong et all 
20041 ; I Yuan. Ouataert & Naravarj|2004K and is also predicted by 



some theoretical model calculations jAbramowicz et alJ fl995 
iLiu et aljfl999l : iRozanska & Czemy||200d : ISpruit & Deufelll2002l) . 

In all these models, the transition radius r tr (r tr = Rti/Rs, and 
Rs = GMbh/c 2 ) is expected to increase with decreasing accretion 
rate m. 

The evolution of central engines in AGNs is mainly gov- 
erned by accretion processes. The black hole accretion processes 
are regulated by the gases supplied near the black hole, i.e., 
sufficient gases supplie d probably lead to a high accretion 
rate, and vice versa {Springel. Di Matteo & Hernquistl 120051: 
Di Ma tteo. Springel & Hernquistl 120051 ; IHopkins. & Hernquistl 



2006t) . As the gases being swallowed by the black hole, the 
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accretion rate decreases with time. An AGN will die off, while 
the gases near the black hole are finally exhausted. How the 
accretion rate m evolves with time remains to be an open issue, 
though an exponentially time-dependent accreti on rate rh(t) was 
widel y employed in ma n y prev io us works (e.g., Park & Vi shniad 
199d: IHaiman & Loebl Il998l : iKauffmann & Haehneltl l2000h . 
McMi llan. Lightman & Cohnl (1981) found accretion rate 



rh(t) oc t~ 2 , if the gases accreted by the black hole are assumed 
to be supplied by stellar collisions or tidal disruptions in a 
dense star system surrounding the black hole. Compared with 
exponentially time-dependent accretion rate, this form of accretion 
rate changes very slowly with time. Many quasar evolution 
model calculations showed that the exponentially time-dependent 
quasar light curve (or simply a step function quasar light curve) 
can well reproduce the obs e rved quasar luminosit y functions 
(e.g., IHaiman & Loebl 1 19981 : IWvithe & Loebl |2002|) . Recently, 
the numerical simulations on the quasar activity triggered by 
the galaxy merger sho wed that the quasar accretion rat e curve 
is very complicated {Springel. Di Matteo & Hernquistl 120051 ; 
|Pi Matteo. Springel & Hernquistl 120051) . In their simulations, the 
gases near the black hole are blown away by the bright quasar 
radiation, and then accretion rate declines rapidly to switch off the 
quasar activity. 

In principle, the evolution of the AGN light curve 
L(t) can be derived from the Eddington ratio distribu- 
tion of an AGN sample, because the Eddington ratio dis- 
tribution £[loK(L/L Edd ) 1 oc [d loff(L/L Ed d)/dtr* (e 
Begel man & Celottil |2004| ; iHopki ns. Naravan & Hernquistl uOO 
However, it is still not straightforward to derive the time-dependent 
accretion rate rh(t), because the radiative efficiency is no longer a 
constant when the accretion disk is evolving from slim disk (high- 
m) (e.g., Abra mowicz et all 1 19881 ; IWang et alj [l999) to a RIAF 



(low-rh) (e.g jNaravan & Yilll994h 



In this work, we derive how the accretion rate m(t) evolves 
with time from an observed Eddington ratio distribution for a low- 
luminosity AGN sample in the local universe based on the stan- 
dard disk/RIAF transition model. In our calculations, the physics 
of different accretion models and the accretion mode transition is 
properly considered. 



2 SPECTRA OF RADIATIVELY INEFFICIENT 
ACCRETION FLOWS 

The transition of a standard thin disk to a RIAF occurs, when the 
accretion rate m decreases to a value below m cr i t . The structure 
of the RIAF is well described b y the self-similar s olu tion, except 
the region near the black hole dNaravan & Yilll995l ; Pvll ll996h . The 
spectrum of a RIAF, L\ (Mbh, rh, a, /3), can be calculated approx- 
imately based on the se l f-similar solution (e.g.. rMahadevarif fl997l ; 
IChang. Choi & Yill2002hlWu & Cadl2005l) . if the parameters M bh , 
m, a, and the fraction of the magnetic pressure /3, are specified. 
However, the physical quantities of a self-similar RIAF deviate sig- 
nificantly from the global solution of the RIAF in its inner region, 
which may lead to inaccuracy for its spectral calculation, because 
most gravitational energy of the accreting matter is released in the 
inner region of the accretion flow. 

In this work, we calculate the global structure of the accretion 
flows surrounding massive Schwarzschild black holes. The accre- 
tion flow is descri bed by a set of g eneral relativistic hydrodynami- 
cal equations (see lManmotcll2000l . for the details). All the radiation 
processes are included in the calculations of the global accretion 



flow structure. Integrating these equations from the outer boundary 
of the flow at R = i? ou t inwards the black hole, we can obtain 
the global structure of the accretion flow passing the sonic point 
smoothly to the black hole horizon. In our calculations, the values 
of som e parameters adopted are different from those in lManmotd 
(2000), which will be discussed in Section 7. We can then calculate 
the spectrum of the accretion flow based on this global structure of 
the RIAF. In the spectral calculations, the gravitational redshift ef- 
fect is considered, while the relativistic optics near the black hole is 
neglected. We only calculate the total luminosity radiated from the 
RIAF without considering its inclination. The derived spectrum in 
this way can be taken as an average spectrum for AGNs, which is 
a good approximation, as AGNs should have randomly distributed 
orientations. 



3 SPECTRA OF STANDARD DISKS 

For a standard t hin disk, the flux due to v iscous dissipation in unit 
surface area is JShakura & Sunvaevll 19731) 



F viB (R) 



3GM bh M 
8-kR 3 



1/2' 



(1) 



where Rs = 2GMbh/c 2 . The local disk temperature of the thin 
cold disk is 



Tdisk(-R) 



1/4 



(2) 



by assuming local blackbody emission. In order to calculate the 
disk spectrum, we include an empirical color correction for the disk 
therm al emission a s a function of radius. The correction has the 
form dChiangll2002l) 



/col (Tdisk) — fcx 



(/« - l)[l+exp(-Vb/Ai/) 



(3) 



1 + exp[(^ p - Vb)/Av] 

where v v = 2.82fc_BTdi s k//i is the peak frequency of blackbody 
emission with temperature Tdj a k. This expression for / co ] goes 
from unity at low temp eratures to fog at high temperatures with 
a transition at vt, ~ u p . IChiand J2OO2I) found that /oo = 2.3 and 
Vb = Av = 5 x 1 15 H z can well reproduce the model disk spectra 
of Hu benv et alj d200 lh . The disk spectra can therefore be calcu- 
lated by 



oc 



rdr 



/col[ eX P(W/colfcsTdisk) - 1] 



,(4) 



where n n = R[ n /Rs is the inner radius of the standard disk. At 
a high accretion rate, m > m cr it> the standard thin disk extends 
to the minimum stable orbit of the black hole, n n = 3, for a 
non-rotating black hole. For a RIAF+standard thin disk system, the 
spectrum emitted from the standard disk region can be calculated 
by using the transition radius r tr (rti = Rti/Rs ) instead of r ln as 
the lower integral limit in Eq. ©■ 



4 TRANSITION RADIUS Rtr BETWEEN RIAF AND 
STANDARD THIN DISK 

In this work, we assume the transition from a standard thin disk to 
a RIAF to o ccur whenever m < m cr j t , i.e., so-calle d "strong prin- 
ciple" (e.g., lNaravan. M ahadeva n~& Ouataertll 19981) . The RIAF is 
naturally expected to match a standard thin disk at the transition 
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radius r tr . The detailed physics, causing such a transition of a stan- 
dard thin disk to a RIAF, is still unclear. It is suggested that the 
standard thin disk transits is truncated at an initial transition radius 
rtr.o an d a RIAF is present within this radius, when m = m cr it- 
The transition radius r tT increases with decreasing accretion rate m 
as 



r tr oc m 



(5) 



where p — 2 is predicted, based on the scenario of transitio n 
triggered by the thermal instability, by I Abramowicz et alj dl995h : 
or p ~ 0.8 — 1.3 i s expected by the disk evaporation induced 
transition scenarios |Liu et al.l fl999: Rozansk a & Czernvl l200d : 
ISpruit & Deufeihrxm In either one of these transition scenarios, 
the transition radius r tr always increases with decreasing accretion 
rate m. The precise initial transition radius rt r ,o is still unknown, 
though it should be small. In this work, we adopt r tr ,o = 20 in all 
our calculations. 



5 THE OBSERVED DISTRIBUTION OF THE 
EDDINGTON RATIO L B /Ledd 

In this work, we adopt the sample given by iHol J2002h . in which 
74 nearby supermassive black holes with both measured masses 
and B-band luminosities. The black hole masses of the sources 
in this sample have been measured by using two different ap- 
proaches: stellar and gas kin ematics and reverberation mapping. 
iTerashima. Ho & Ptakl j2000l) found that X-ray luminosities in 2- 
10 keV of LINERs (low-ionization nuclear emission-line regions) 
with broad Ha emission in their optical spectra are proportional 
to their Ha luminosities. This indicates that the dominant ioniz- 
ing source in LINERs is photoionization by hard photons from 
low-luminosity AGNs. The B-band luminosities of the sources in 
this sample are estimated from the line emissions, which are sup- 
posed to be photo-ionized by the nuclear radiations. This may cause 
some uncertainties for individual sources, but the derived Edding- 
ton ratio should still be reliable in statistical sense. This sample 
includes 17 PG quasars, which may not be in the same popula- 
tion with those low-luminosity counterparts and are accreting at 
higher rates m > m C rit. As we are focusing on the inefficient ac- 
cretion history of AGNs, we leave out these 17 PG quasars, which 
leads to 57 sources. The Eddington ratio distribution for this sam- 
ple is plotted in Fig. [T] This samp le was used to explore the evo- 
lution of low-luminosity AGNs by Hopkin s. Naravan & Hernquist 
(2006). The selection effects of this sample were extensively ana- 
lyzed in their work (see Section 2.1 in their paper for the details). 
They found that, the FR I radio galxies from the sample used by 
iMarchesini. Celotti & Ferrareseet fc004l) exhibit significant differ- 
ent luminosity distribution from that of the sample bv lHol j2002t) . 
while both of these two sample exhibit similar Eddington ratio dis- 
tributions. This means that the results of AGN evolution derived 
from either one of these two samples would be qualitatively un- 
changed, which may imply that the derived results are not affected 
significantly by selection effects. 

The maximal distance for the sources measured by reverbera- 
tion mapping method is about 3 ti mes that of the subsample mea- 
sured by the stellar/gas kinematics. Hopki ns. Naravan & Hernquist 
(2006) suggested to multiply the relative fraction of the sources 
measured kinematics by 3 3 to account for the different subsam- 
ple volumes. The derived distribution for this tentatively volume- 
corrected sample is also plotted in Fig.Q] 



0.3 
0.25 



° 0.2 

o 
CO 

E 



0.05 



-10 



-8 



-6 



-4 



-2 



'°9io V L Edd 



Figure 1. The distribution of the Eddington ratio L B / ^Edd at B-band for 
the sample given bvlHo1 j2002h (the solid line). The dotted line is the dis- 
tribution for the effective volume-corrected sample (see Section 5 for the 
details). 



6 THE TIME-DEPENDENT ACCRETION RATE M{T) 
DERIVED FROM THE EDDINGTON RATIO 
DISTRIBUTION 

The light curve of AGNs at B-band can be derived from the ob- 
served Eddington ratio distribution diV/dlog(I/ B /LEdd): 



dlog(L B /L E dd) 
At 



TV 1, 



AN 



dlog(L B /L B dd) 



(6) 



where N tot is the total number of the sources in the sample, and t 
is normalized to unity. The time-dependent accretion rate m(t) can 
be derived from the light curve by 



dm(t) _ dlog(L B /L B dd) ( dlog(L B /I/Edd) 



At 



At 



Am 



(7) 



where d log(L B / L-Edd)/Am is available based on spectral calcula- 
tions for the accretion mode transition model desribed in Sections 
2-4. Here, we have to assume that all sources have the same time- 
dependent accretion rate rh(t), and it evolves monotonically with 
time. The latter assumption may not be the case for some individual 
sources in short timescales, but it should be reseasonable in statistic 
sense for a sample of AGNs, because accretion rates should decline 
in a long timescale. Substituting Eq. $6$ into Eq. I0, we obtain 



(8) 



Am(t) _ Ntot (ANy 1 
dt \Am) 

Using the X-ray luminosity function of AGNs given bv lUeda et al.l 
(2003), we calculate number density of AGNs in comoving space 
as a function of redshift. It is found that the number density of 
AGNs in the local universe is about 45 per cent of that at z = 0.5, 
which implies that the AGNs at low redshifts are roughly in a steady 
evolving state with switching on and off being in balance. The cos- 
mological evolution of AGNs at low redshifts should be unimpor- 
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in 



Figure 2. The spectral evolution of RIAF+standard thin disk systems with 
varying transition radii rt r . The dotted lines represent the emission from 
the outer standard thin disk regions. The solid line represents the case for 
p = 1, while the dashed line is for p = 2. 

tant, which will affect our results very little, as the sample used in 
this work is limited to the sources in the local universe. 



7 RESULTS 

The detailed physics of the viscosity in accretion disks is still 
quite unclear, and it is usually described by the viscosity param- 
eter a. Assuming accretion to be driven by turbulent stresses gen- 
erated by the magnetorotational instability, the three-dimensional 
MHD simulations suggest t hat the viscosit y parameter a in the 
accretion flows is ~ .1 jArmitagel 1 1998b . or ~ 0.05 — 0.2 
dHawlev & Balbuj I2002T) . The critical accretion rate rh cr it ~ 
0.01 for accretion mode transition is suggested by different au- 
thors either from observations or theor etical arguements (see, e.g., 
Naravan, Mahad evan & OuataerJ Il998h . Our numerical calcula- 
tions for the global structure of the flows show m cr i t — 0.01 
for q = 0.2. In this work, we adopt the a-viscosity, a = 
0.2, and limit our calculations of the RIAF structure to hydro- 
dynamics, but the magnetic pressure is included by a parameter 
f3. The parameter f), defined as p m = (1 — 0)pto%, describes 
the magnetic field strength of the gases in the flow. This pa- 
rameter is in fact not a free parameter, which is related to the 
viscosity parameter a as (3 ~ (6a — 3)/(4q — 3), as sug- 
gested by the MHD simulations lHawlev. Gammie & Balbusl 19961 ; 
iNaravan. Mahadevan & Ouataertlll998l) . For a = 0.2. Q ~ 0.8 is 
required. iManmotol (2000) adopted very small values of 8 in the 
calculations, as those of tra ditional advecti on dominated accretion 
flow (ADAF) models (see iNaravarJ |2002| . for a review, and ref- 
erences therein). It was pointed out that such a small 8 adopted 
in traditional ADAF models is very unlikely, because a signifi- 
cant fraction of the viscously dissipated energy (up to 8 ~ 0.5) 
could go into electrons by magn etic reconnection, if the magnetic 
fields in the flow are strong dBisnovatvi-Kogan & Lovelace! 1 19971 
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Figure 3. The derived time-dependent accret ion ra tes rh(t). The blue lines 
represent the cases derived directly from Ho l2002h 's sample, while the red 
lines are for the effective volume-corrected sample. The solid lines are for 
the cases of p = 1, while the dotted lines are forp = 2. 

2000). The average black hole mass for the sources in this sam- 
ple is ~ 2.4 x 10 8 Mq. In our calculations, we adopt 8 — 0.5, 
and Mhh = 2.4 x 1O 8 M0. Our numerical calculations show that 
the dimensionless light curves Lb /L-^dd as function of rh change 
very little for different black hole masses adopted (in the range of 
~ 10 6 ~ 10 M Q ). 

Using the method described in Sections 2 and 3, we can cal- 
culate the global accretion flow structure and then the spectrum 
varying with accretion rate rh and r tr . The initial transition radius 
rti : o = 20 is adopted in all our calculations for RIAFs. The Ed- 
dington ratios Lb / LEdd evolving with m are plotted in Fig. [2] The 
time-dependent accretion rates rh(t) can be derived from the ob- 
served Eddington ratio distributions by using Eqs. l[6} and ([7J (see 
Fig.0. 



8 DISCUSSION 

For a RIAF+standard thin disk system, its optical emission is 
mostly from the outer thin accretion disk region, if its accretion rate 
rh is close to the critical value merit, because the transition radius 
is small for this case. For a large transition radius, corresponding 
to a low rh, the temperature of the standard thin disk region is very 
low, and the optical emission is dominantly from the inner RIAF 
(see Fig.|2]l- The optical luminosity drops rapidly with rh, when the 
flow is accreting at < rh cr it- 

In Fig. [2] we find that the luminosity decreases rapidly after 
the accretion mode transition at rh — rh cr it ~ 0.01. It is found that 
the light curves LB(rh) are almost same for p = 1 and p = 2, ex- 
cept the accretion rate rh is in the range of ~ 0.004 — 0.01. When 
the accretion rate rh is slightly lower than rh cr i t , the emission is 
dominantly from the outer thin disk regions, because the standard 
thin disks still extend to small radii. The emission becomes dom- 
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inated by the radiation from the inner RIAFs, when the accretion 
rate m < 5 x 10" 3 . 

Based on the relations between accretion rate rh and B- 
band luminosity Lb (rh) given in Fig. [2] we can derive the time- 
dependent accretion rates rh(t) from the observed Eddington ratio 
distributions for AGNs. In Fig. [3] we plot the derived rh(t) for the 
cases with two different values of p for variable transitio n radi us 
models, from the Eddington ratio distribution of the IHqI 
sample or the sample with effective volume-corrected (see Section 
5), respectively. We find that there is a rapid decrease of rh(t) from 
10" 2 to ~ 10~ 4 either for the cases with different values of p or 
from different Eddington ratio distributions (the original one or the 
effective volume-corrected one). It is not surprised to find that the 
derived time-dependent accretion rates rh(t) are quite similar for 
the cases with p — 1 or p = 2, because the theoretical curves 
Lb(»ti) are very similar for these two different values of p (see 
Fig-O- F° r me Eddington ratio distribution of the effective volume- 
corrected sample, there are more sources with low luminosities, i.e., 
less sources with rh close to rh cr i t , so the derived time-dependent 
rh(t) decreases more rapidly than that derived from the original 
sample. We find that the main feature of the time-dependent rh(t) 
has not bee n changed by this effective volume-corrected sample. 
IcHol l l2005t) calculated the hard X-ray emission from all RIAFs in 
faint AGNs, and compared them with the observed X-ray back- 
ground. It was found that the accretion rate should decrease rapidly 
to the value far below the critical rate within a timescale shorter 
than 5 per cent of bright quasars' lifetime, which is consistent with 
our present time-dependent accretion rates rh(t) derived from the 
Eddington ratio distributions. The present derived time-dependent 
rh(t) is based on the assumption of monotonical evolution with 
time. From Eq. l[8}, we find that drh(t) /At also represents the dis- 
tribution of accretion rates rh for the sample. 

The present sample may have overlooked a number of sources 
with very low Eddington ratios (< 10 -5 — 10 -6 ), which means a 
complete sample should include more sources accreting at very low 
rates < 10~ 4 than the present sample. This implies that the relative 
timescale of the rapid drop of accretion rate from 0.01 to ~ 10~ 4 
should be even shorter than those plotted in Fig. [3] if a complete 
sample is adopted. Such a complete sample is still unavailable now, 
however, the main feature of rh(t) with a rapid declining between 
0.01 and ~ 10 -4 derived in this work will not be changed qualita- 
tively. 

The fraction of the viscously dissipated energy that directly 
go es into electrons, S = 0.5, is adopted i n this work, as suggested 
by iBisnovatvi-Kogan & Lovelace! dl997t) . The precise value of S 
is still unknown, and it may slightly be lower than 0.5. If a slightly 
lower 5 is adopted, the derived light curves Lb (rh) have the similar 
form, but they should be systematically lower than those in Fig. [2] 
Thus, the main feature of rh(t), a rapid decrease after the accretion 
mode transition, will still be present as those plotted in Fig. [3] 

The RIAF may have winds, and a power-law r-dependent ac- 
cretion rate is assumed, thou gh the detailed physics is still unclear 
jBlandford & BegelmarJI 19991) . When the accretion rate rh is close 
to 0.01, the optical emission is dominated by that from the outer 
standard thin disk region (see Fig. O. Thus, the optical spectrum 
of the RIAF+standard thin disk system will not be affected by the 
winds of the RIAF, when rh is close to 0.01. The optical emission 
from the RIAFs is attributed to Compton up-scatterings of the soft 
synchrotron photons by the hot electrons in the accretion flows. Our 
numerical calculations on the RIAF spectra show that most optical 
emission (more than 90 per cent) is from the regions within 10 Rs- 
This implies that the optical emission of a RIAF is similar to a 



RIAF with winds, provided they have the same accretion rate at 
their inner edge. In this work, we focus on the accretion rate at the 
inner edge of the accretion flow, i.e., the mass of gases swallowed 
by the black hole. Thus, the derived time-dependent accretion rates 
at the inner edges of the RIAFs will not be altered much, even if 
winds are present in the RIAFs. 

In all our calculations, the initial transition radius r tr ,o = 20 is 
adopted. As the transition radius r tr increases rapidly with decreas- 
ing accretion rate rh (see Eq. [5J, the resulted light curves Lb (rh) 
will only be slightly different for different values of r tr ,o adopted, 
provided r tr ,o is not very large. Thus, the derived rh(t) is insensi- 
tive to the value of rtr.o. 

In principle, one has to consider the activity trigger rates of 
galaxies along the cosmic time in the study of AGN evolution. 
However, it is difficult to have a complete sample including both 
bright and faint AGNs, which can be used to explore the evolution 
of bright quasars to faint AGNs. In this work, we only explore the 
time-dependent rh(t) for rh < m cr i t . and the sources in the sample 
we used are limited in the local universe. Thus, the results may not 
be affected by the unknown trigger rates, though it is unclear if the 
derived rh(t) in this work is valid for the AGNs at high redshifts. 
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